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ABSTRACT 

A soft ultraviolet radiation field, 10.2 eV< hv < 13.6 eV, that permeates neutral 
intergalactic gas during the Epoch of Reionization (EoR) excites the 2p (directly) 
and 2s (indirectly) states of atomic hydrogen. Because the 2s state is metastable, 
the lifetime of atoms in this level is relatively long, which may cause the 2s state 
to be overpopulated relative to the 2p state. It has recently been proposed that for 
this reason, neutral intergalactic atomic hydrogen gas may be detected in absorption 
in its 3-cm fine-structure line (2siy 2 - ► 2p 3 / 2 ) against the Cosmic Microwave Back- 
ground out to very high redshifts. In particular, the optical depth in the fine-structure 
line through neutral intergalactic gas surrounding bright quasars during the EoR may 
reach Tps ~ 10~ 5 . The resulting surface brightness temperature of tens of /i K (in ab- 
sorption) may be detectable with existing radio telescopes. Motivated by this exciting 
proposal, we perform a detailed analysis of the transfer of Ly/3, 7, 5, ... radiation, and 
re-analyze the detectability of the fine-structure line in neutral intergalactic gas sur- 
rounding high-redshift quasars. We find that proper radiative transfer modeling causes 
the fine-structure absorption signature to be reduced tremendously to Tpg <10 -10 . 
We therefore conclude that neutral intergalactic gas during the EoR cannot reveal its 
presence in the 3-cm fine-structure line to existing radio telescopes. 
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1 INTRODUCTION 

It has recently been proposed that it may be possible to 
observe neutral atomic hydrogen gas during the Epoch of 
Reionization (EoR), in absorption against the Cosmic Mi- 
crowave Backgr ound (CMB) in its 3-cm fine-structure line 
(2s 1/2 -> 2pq/,. ISethi et al.ll2007l . hereafter S07). 

The possible detectability of the 3-cm line derives from 
the fact that the 2s s tate of atomic hydrogen is metastable 
|Breit fc Tellerlll94d ). That is, the average lifetime of an 
undisturbed atom in the 2s-state is tzs = l/^bsis ~ 0.1 
s, aft er which it decays into the Is state by emitting 2 pho- 
tons l|Spitzer fc Greenstein|[l95l! ) . For comparison, the mean 
life-time of atoms in the 2p state is t2 P = 1 /A^pis ~ 10 -9 s. 
Hence, when atoms are excited into their 2s and 2p states at 
comparable rates, the 2s level is expected to be highly over- 
populated relative to the 2p level. This may result in a non- 
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negligible op tical depth in the 2s!/ 2 — 2p 3 / 2 fine-structure 
line (see e.g. IWild|[l952l) . 

The detectability of this line has been discussed in 
the context of HII regions in the local Universe (e.g. 
iDennison et alj|2005l . for a more detailed discussion). In lo- 
cal HII regions, the 2s level is populated mainly following 
recombination, while both recombination and photoexcita- 
tion by Lya photons populate the 2p state. The Lya excita- 
tion rate of the 2p-state is limited by dust, and the 2s-state 
is typically heavily overpopulated in HII regions which re- 
sult in 3-cm optical depths of order tfs ~ 10 -3 (Ershov 
1987; Dennison et al, 2005, where the absorption is mea- 
sured relative to the free-free continuum that is produced 
in the HII region). S07 have argued that during the EoR, 
bright quasars that are luminous in the restframe soft-UV 
(10.2 eV < < 13.6 eV) may indirectly photoexcite the 2s 
state of neutral hydrogen gas in the low density intergalactic 
medium to levels such that tfs ~ 10 -5 for CMB photons 
passing through these region^. Absorption of these CMB 
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Figure 1. Schematic diagram of the energy levels of a hydrogen 
atom. On the left side, we show the first 3 energy levels, and 
their splitting into levels with different angular momentum. The 
notation for each level is nL, where n is the principle quantum 
number, and L denotes the orbital angular momentum number 
(L = S corresponds to I = 0, L = P corresponds to I = 1, L = D 
corresponds to I = 2). The dotted lines show the allowed transi- 
tions through the emission or absorption of a single photon. In 
a 3-level atom, the selection rules only allow the 2s state to be 
populated by absorbing a Ly/3 photon, and by subsequently emit- 
ting an Ha photon. On the right, the fine+hyperfine-structure 
splitting of the 2s and 2p levels is shown in more detail. The 
Lamb-Rutherford shift causes the 2s 



1/2 



level to lie above the 



2p 1/2 -lcvcl by AE/hp ~ 1.1 Ghz (e.g. 1 Wild|[l953 ) . The 3 allowed 
transitions between the 2si/2 and 2p3/2 levels are indicated with 
dashed lines, and the frequency (in M Hz) of each transition is 
indicated (also see lDennison et alj|2005l ). The relativ e strength of 
each transition 9852:9876:10030 is 1:5:2 d Wildlll952h . 



photons would result in an absorption feature with a bright- 
ness temperature of a few tens of pK at 1.4 Ghz, which may 
be detectable with existing dish antennas such as the Green 
Bank Telescope (GBT), Parkes and the most compact con- 
figuration ofthe Australia Telescope Compact Array (ATCA, 
e.g.lCarilli 2008:). This result is important, because it implies 
that neutral intergalactic atomic hydrogen gas that exists 
during EoR may reveal itself in a line that is different (and 
possibly easier to d etect) than the well s tudied 21-cm hyper- 
fine transition (e.g. IScott fc Rees|[l990l ; lLoeb fc Zaldarriagal 



12004 iFurlanetto et al Jl2006t l . Motivated by this exciting re- 
sult, and by observational efforts to detect this transition, 
we re-analyze the detectability of the fine-structure line dur- 
ing the EoR. We perform more detailed modeling of the 
transfer of Ly/3 radiation than originally discussed by S07. 
We will show that (unfortunately) proper modeling of Ly/3 
(and Ly7, 6, ...) radiative transfer radically lowers the value 
of tfs to levels that cannot be probed with existing radio 
telescopes. 

The outline of this paper is as follows. We review the 
dominant excitation mechanisms of the 2s and 2p levels in 
Ej2j In |3l we compute the effective excitation rate of the 
2s-state during the EoR, by properly accounting for ra- 



of atomic hydrogen can produce a detectable opacity. One can 
estimate Tps by comparing it to the Gunn-Peterson optical depth 
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10 5 for a fully neutral IGM at z = 6, rp$ = 10 -5 only requires 
that ~ 10 — 12 . It will be shown in § l5,ll that these levels can 
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be reached in regions around high-redshift quasars. 



diative transfer, and discuss the detectability of the 3-cm 
fine-structure line in [jH The validity of our model assump- 
tions, and the implications of this work are discussed in 
Sj5] Finally, we present the conclusions in |6] The param- 
eters for the background cosmology used throughout are 
(Q m n A> n t ,fe) = (0.27,0.73,0.042,0.70) (jKomatsu et all 
12008}) . 



2 THE EXCITATION MECHANISM OF THE 
2S LEVEL 

Figure [1] depicts the first three principle energy levels of 
the hydrogen atom. On the left, fine and hyperfine-structure 
splitting of the lines are ignored. The notation that is used is 
nL, where n is the principle quantum number, and L denotes 
the orbital angular momentum number (L — s corresponds 
to I — 0, L — p corresponds to I = 1, L = d corresponds to 
I = 2). In most astrophysical conditions, the vast majority 
of hydrogen atoms are in the ls-state. Atoms can be excited 
out of the ground state through the following processes: (i) 
collisional excitation, (n) recombination following photoion- 
ization or collisional ionization, and (Hi) photoexcitation. 
Following S07, processes (i-ii) are ignored (see § 15.11 - under 
roman numeral I -for a more detailed motivation). 

The quantum mechanical selection rules state that the 
only atomic transitions that are allowed by emitting or ab- 
sorbing a photon are those in which I changes by Al = ±1. 
Hence, absorption only results in transitions of the form 
Is — > np, while subsequent emission only results in np — > 
ms. The allowed transitions among the first 3 levels of the 
H-atom are indicated in Figure [T] as dotted lines. In a 3-level 
atom (the impact of using a multilevel atom is discussed in 
§ I5.1l under roman numeral II) the selection rules only allow 
the 2s state to be populated by absorbing a Ly/3 photon, 
and by subsequently emitting an Hcv photon. Hence, the 
level population of atoms in the 2s state is determined by 
the effective Lyf3 scattering rate, denoted by Pp, e s, which is 
defined as the rate at which Ly/3 photons absorbed, and sub- 
sequently converted into Ha + 2 continuum photons. This 
rate is lower than to the "ordinary" Ly/3 scattering rate, Pp, 
which is the rate at which Ly/3 photons are absorbed. 



3 THE 2S AND 2P LEVEL POPULATIONS 

As mentioned above, the effective Ly/3 scattering rate, P/3, e ff , 
is crucial in calculating the number of atoms in the 2s-state. 
The goal of this section is to evaluate P/3, e ff(r). In § 13.1 
we describe the Monte-Carlo technique that was used to 
compute the general and effective Ly/3 scattering rates, and 
in § 13.21 we present our results. 



3.1 Outline of Monte-Carlo Calculation 

We compute the general and effective Ly/3 scattering rate 
by performing Monte-Carlo (MC) calculation of the Ly/3 ra- 
diative transfer in the neutral gas surrounding a spherical 
HII region of radius i?Hii = 5 Mpc, that is based on the 
physical size of the HII region that possibly exist around 
z = 6 quasars l|Mesinger fc Haimanl [20o3 : IWvithe fc Loebl 
12004 but see Lidz et al 2007, Bolton & Haehnelt 2007). 
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This MC-approach to radiative transfer, in which the tra- 
jectories of individ ual photons are followed, is common for 
Lyg radiation (e.g. ILoeb fe Rybickil Il999l : lAhn et alJbood : 
I Zheng fc Miralda-Escudel 120021). We use th e Lya-MC code 
originally described in Diikstra et al, (2006), which is easily 
adapted for Ly/3 radiative transfer. The calculation is per- 
formed as follows: 

1. A photon is emitted with a frequency that is drawn 
randomly from the range u £ [0.95^, 1.05^]. The transfer 
calculation therefore focuses on photons in a frequency range 
Ai/p — Q.lvp centered on up, similar to the optically thin 
calculation that was originally performed by S07 (and which 
is discussed in Appendix IA1I) . 

2. The photon propagates an optical depth r = — log R, 
where R is a random number between and 1. The 
optical depth translates to a physical distance via r = 
f^ds n m (s)ap(u[l - Hs/c]). 

3. Each time a Ly/3 photon is absorbed, there is a 
P/3,scat = 0.88 probability for the atom to decay directly 
back to the ground state by emitti ng a Ly/3 photon (jHiratal 
120061 ; iPritchard fc FurlanettoboOot ) . This is simulated in the 
MC-calculation by generating a random number R £ [0, 1]. 
The photon scatters as a Ly/3 photon if R ^ p/3, SC at, in which 
case a new direction is chosen and we return to step 2. On 
the other hand, a transition of the form 3p — > 2s occurs when 
R > P/3, scat- After this process, the Ly/3 photon is destroyed, 
and we return to step 1.. 

We compute total Ly/3 scattering rate by keeping track 
of the number of scattering events that occur in shell j, 
Nj (in the MC-simulation, the gas is sampled by a dis- 
crete number of gas shells, see Dijkstra et al 2006). The 
total Ly/3 scattering rate transfer in shell j follows from 
Pp,j = x Nj/(N-H,j x Nm C ), where Np is the total rate 
at which the quasar is emitting photons in the frequency 
range v £ [0.95^,3, 1.05i//j] (we have assumed Np = 10 58 s _1 
to compare to S07 as closely as possible), AThj is the to- 
tal number of hydrogen atoms in shell j, and where Af mc 
is the total number of photons that is used in the MC- 
simulation. The effective Ly/3 scattering rate is obtained 
similarly from Pp,j, e s = Np x A r j,3 P 2s/(A i 'Hj x A^ic), where 
Nj t 3 P 2s is the total number of scattering events of the se- 
quence Is 3p 5? 2s in shell j. 



3.2 The General and Effective Ly/3 Scattering 
Rates derived from the MC-Simulation 

The general and effective Ly/3 scattering rates we obtain 
from the Monte-Carlo calculation are shown in Figure [2] as 
the black and red solid histograms, respectively. The blue 
dotted line shows the scattering rate that is obtained when 
the gas is assumed to be optically thin (which corresponds 
to the original calculation of S07, see Appendix l Ali tor more 
details). Furthermore, the blue dashed line shows the scat- 
tering rate that one obtains if one ignores that Ly/3 photons 
undergo multiple scattering events, and simply suppresses 
the incoming Ly/3 flux by a factor of e~ T (this calculation is 
described in Appendix I A2|) . 
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2 The scattering process is assumed to be partially coherent, and 
described by a Rayleigh phase function, (see Dijkstra & Loeb 2008 
for a more detailed discussion of the scattering phase function). 



Figure 2. The effective (general) Ly/3 scattering rate that was 
derived from the MC-simulation is plotted as the red (black) his- 
togram. The inset zooms in on the edge of the HII region. Also 
shown as the blue dashed line is the scattering rate that one ob- 
tains by suppressing the Ly/3 mean intensity by a factor of e — T 
(§E2J. The green dot-dashed line (inset only) includes the contri- 
bution from higher Lyman series lines (see jjSTT) . The dotted line 
shows the scattering rate in the optically-thin limit. This Fig- 
ure shows that: (4) the effective Ly/3 scattering rate derived from 
the MC-simulation is reproduced well by simply suppressing the 
mean Ly/3 intensity by a factor of e~ T (see § IA2I I: and (ii) the 
effective scattering rate decreases rapidly with radius, (see inset) 
and levels off at a level that is ~ 5 orders of magnitude smaller 
than the rate obtained from the optically thin calculation. Lastly, 
the green dotted line in the main plot shows the scattering rate 
of photons that would have redshifted into resonance at radius=8 
Mpc (see text) 



Figure [5] shows that the effective Ly/3 scattering rate 
is much lower than the rate that is obtained in the opti- 
cally thin calculation. At small distances into the HI region 
(d > 10 -2 kpc, which is shown in the inset), the effective 
scattering rate is Pp, c a < 10~ 13 s _1 . The effective scatter- 
ing rate decreases rapidly with radius, and levels off at a 
level that is ~ 5 orders of magnitude smaller than the rate 
one obtains from the optically thin calculation. This rapid 
decrease owes to the large line-center optical depth of the 
neutral IGM to Ly/3 photons, which causes resonant scat- 
tering only to occur near the edge of the HII region. This 
point is illustrated by the green dotted line in the main plot 
of Figure [21 which shows the scattering rate of photons with 
frequencies such that they would have redshifted into res- 
onance at r = 8 Mpc (i.e. 3 Mpc into the neutral IGM). 
When these photons reach r — 7 Mpc, they find themselves 
~ 750 km s blueward of the resonance, and their absorp- 
tion cross-section is suppressed by a factor of ~ 10 -8 rela- 
tive to the resonance absorption cross-section. The scatter- 
ing rate of these photons is therefore also suppressed by a 
factor of 10 -8 (relative to the optically thin rates). The scat- 
tering rate increases rapidly with radius as the photons red- 
shift closer to resonance, and the maximum scattering rate 
occurs at r = 7.95 Mpc, where the photons have redshifted 
to within ~ 3 Doppler widths from resonance. The optical 
depth of the IGM increases rapidly as these photons red- 
shift closer to resonance, which results in a sharp cut-off in 
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the scattering rate at r >7.95 Mpc. In other words, photons 
that were emitted somewhat blueward of the Ly/3 resonance 
frequency -which need to travel far into neutral gas from 
the HII region edge before redshifting into the line center- 
have a high probability of being absorbed in the damping 
wing of the line, and do not generally make it to line cen- 
ter. The absorption cross section in the damping wing is low 
compared to that at line center, and so these photons yield 
a correspondingly small scattering rate. Therefore the Lyj3 
scattering rates are strongly suppressed at d J>10 -2 kpqj. 
The sharp cut-off of the effective Ly/3 scattering rate is also 
discussed in Appendix I A2I 

Furthermore, Figure [2] shows that the effective scatter- 
ing rate that is extracted from the MC-simulation (red solid 
histogram) is well reproduced by the model in which the 
mean Ly/3 intensity is suppressed by a factor of e _T (the blue 
dashed line). The reason for this is that each Ly/3 photon can 
only excite 1 atom into the 2s state: each time a transition of 

the from Is L -^+ 3p ^? 2s occurs, the Ly/3 photon is lost per- 
manently. This almost corresponds to a calculation in which 
one ignores scattering and simply suppresses the intensity 
with a factor e _T . The only difference is that Ly/3 photons 
typically scatter (N sca .t) = 1/(1 — P/3,scat) ~ 8.3 times be- 
fore being converted into Ha + 2 continuum photons. The 
random walk associated with the scattering process trans- 
ports the Ly/3 photons on average ~ *J (-/V scat ) ~ 2.9 mean 
free paths through the neutral IGM, before indirectly excit- 
ing an atom into its 2s-state. In practise, this introduces a 
difference with the scattering rate obtained from the e~ T - 
calculation only at the smallest separations from the edge 
of the HII region (r < 10 -2 kpc, see inset). At these small- 
est radii, the gas is optically thin to Ly/3 photons, and here 
the total Ly/3 scattering rate (black solid histogram) as ex- 
tracted from the MC-simulation agrees best with the an- 
alytic solution given by Eq (IA6() with r = 0. Naturally, 
here the effective scattering rate is suppressed by a factor of 

(1 - P/3,scat) = 0.12. 



4 THE SIGNAL 

The change in brightness temperature for CMB radiation 
passing through the neutral gas with an enhanced popula- 
tion of atoms in the 2s-state is given by 

r ~r z 

where 7cmb = 2.73(1 + z) is the CMB temperature at red- 
shift z, and T cx is the excitation temperature of the 2s- 
2p transition, which is defined as T ox = \ FS S/2kB, where 

5 — 2hv FS n u j g standard source function (see 

iRvbicki fc Lightmanl 1 19791 . Eq 1.79), ni (n.) denotes the 

3 Note that in our model the total scattering rate, which is given 
by Vtot = nu Ir n dr 4irr 2 Pp tC s( r ) = Np> i- e - photons are con- 
served. On the other hand, Vtot 3> N[j in the optically thin limit 
(JSU- Each emitted Ly/3 photon can only undergo one effective 
scattering event, and the total effective scattering rate can never 
exceed the total rate at which Ly/3 photons are emitted. The op- 
tically thin calculation is therefore logically inconsistent. 



number density of atoms in the 2s level (2p level), and gi (g u ) 
denotes the statistical weight of this level. Throughout we ig- 
nore that hyperfine splitting of the 2s and 2p levels actually 
results in 3 possible 2s-2p transitions (see Fig [1} . We dis- 
cuss the impact of hyperfine splitting on our results in § 15.11 
(under roman numeral V). Furthermo re, the optical depth 
to 2s - 2p absorption, tfs is given by (|Rvbicki fc Lightmanl 
1 19791 . Eq 1.78) 

tfs{v) = J ds a(u(s)) = ^J^^ J ds {^ni~n^j(f)(u{s)), 

(2) 

where Avps = v t h^Fs/ c , B u i is the Einstein B-coemcient 
for the fine-stru cture transition, and 4>(v) is the line pro- 
file function (see lRvbicki fc Lightmanl fl979). S07 move the 
factor — m — n u outside of the integral and arrive at their 

91 ° 

Eq (7). However this is not allowed, because we found Tp(r), 
and therefore n;(= n^s), to be a rapidly changing function of 
radius (see Fig[2}. This makes the evaluation of tfs(j/ ) some- 
what trickier. Fortunately, we may simplify Eq[T]and Eq[2] 
by assuming that ni > n u (which is verified in § 15.11 under 
roman numeral IV). Under this assumption, T ex <C Tomb 
and AT b (u) = — Tfs(i/)Tcmb/(1 + z), and the optical depth 
to 2s- 2p absorption at frequency v simplifies to 

tfs{v) = J ds n2s(s)o-Fs(v[l — Hs/c]), (3) 

where ofs{v) = 4^^"' ^ < t>{ 1 ^)- F° r simplicity, we assume 
that the edge of the HII region is infinitely sharp (i.e. the HII 
region contains no neutral hydrogen, and the neutral frac- 
tion jumps from to 1 instantly at the edge of the HII re- 
gion) and extends out to r = Run = 5 Mpc (see § 15.11 roman 
numeral III). Under this assumption, the optical depth along 
a path that is pointing radially outward from the quasar is 
given by 

/>OC 

tfs{v) = -j^- / dr Pp{r)cr FS {v\\ - Hr/c]), (4) 

A2sls jR mi 

where we used that U2 S = ^ 3 P ^ = ^[f~p ■ The last approxi- 
mation is allowed because the vast majority of all hydrogen 
atoms are in the Is state. 

The life-time of the 2p state is set by the Einstein co- 
efficient of the 2p — > Is transition, A2 P i s = 6.25 x 10 8 
s _1 , and the line profile of the 2s- 2p transition has a Voigt 
parameter a = A 2v i a /(ATrAu FS ) = 117(T gas /10 4 K)~ 1/2 
jPennison et alj 120051 ). The line profile is therefore com- 
pletely dominated by the damping wing, and the cross- 
section displays a much weaker frequency dependence than 
is common among the Lyman-series transitions: e.g. even 
for xfs = 100, where xfs(v) = (y — - v F s) j Av F s , the cross- 
section is only reduced by a factor of <j>(xFs)/4>{fy — 0.57. 

We compute tfs(v) by using Tp(r) we found in § 13.21 
(and we account for the additional contribution of higher 
Lyman series following i]5.1[) into Eq [4] The result of this 
calculation is shown in Figure [3] which shows that the max- 
imum optical depth is tfs, max = 4 x 10~ 10 . This is 5 or- 
ders of magnitude smaller than the value that was found by 
S07. The reason for this large difference is the assumption 

that the factor — rii — n u could be taken outside the inte- 
rn 91 

gral in Eq[2] This is allowed only when the dependence of 
this factor on r is weaker than that of ctfs(i / [1 — Hr/c]). 
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Figure 3. This Figure shows the optical depth in the 2s-2p fine- 
structure transition, Tps(v), through neutral gas surrounding a 
quasar HII region at z = 6.5. The maximum optical depth is 
Tps max ~ 4 X 10 — 10 , which results in a maximum brightness 
temperature of ~ 10 — 10 K. This brightness temperature if five 
orders of magnitude lower than previous estimates, and likely 
precludes detecting hydrogen fine-structure absorption during the 
EoR. 

However, because the 2s- 2p absorption cross-section is com- 
pletely dominated by the damping wing, it has a tremendous 
FWHM that is v FWH m/vfs ~ 0.01. At z = 6.5, Hubble ex- 
pansion over 4 Mpc introduces a comparable frequency shift 
of Av/u ~ 0.01, and S07 therefore implicitly assumed that 
their calculated level population in the 2s-state applied to a 
'shell' that was at least 4 Mpc thick. In sharp contrast, we 
found that their calculation applied only to a shell of gas 
that was ~ 10 -2 kpc thick (see the inset of Fig [2}. In other 
words, the total column of neutral H atoms in the 2s state 
is lower by ~5 orders of magnitude. 

It is worth emphasizing that our calculation of TFS.max 
assumed that the quasar emitted Np = 10 58 photons s _1 
in the frequency range v £ [0.95, 1.05]vp. This is very likely 
a factor of ~ 10 too high (see Appendix IA1|) . Therefore, 
in reality the value of TFS.max is likely even lower by an 
additional factor of ~ 10. 



5 DISCUSSION 

5.1 Validity of Model Assumptions 

(I) : In 52] we mentioned that the collisional excitation 
rate of -and the recombination rate into- the 2s state, were 
negligibly small compared to the effective Ly/3 scattering 
rate. For example, the collisional excitation rate is given by 
C 2s = 2.24xl0- 6 n e r g ~ s 1/2 exp(-10.2 eV/fcT gaa ) (|Osterbrockl 
1 19891 . p54-55). If we assume that the gas is half ionized (n e = 
0.5n H ), then we find that C 2a < 10" 14 s" 1 for T gas < 3 x 10 4 
K. This rate is orders of magnitude less than the effective 
Ly/3 scattering rates in the regions of interest. Similarly, we 
find that the recombination rate is R 2a < Q r cc,Bn e = 10 -17 
s _1 , which is even smaller. 

(II) : So far we have ignored the fact that the 2s state 



of atomic hydrogen can also be accessed following the ab- 
sorption of higher order Lyman series photons, e.g. via a 
scattering event of the form Is 5p — ? 2s. In principle, 
it is straightforward to expand the MC-approach to include 
scattering of these higher-order Lyman series photons, but 
this is not necessary. We found in § 13.21 that the effective 
Ly/3 scattering rate was calculated reasonably well if one 
simply suppresses the incoming Ly/3 flux with a factor of 
e _T . For simplicity, we also use this approach to calculate 
the effective scattering rates of the higher order Lyman se- 
ries photons, and cast the total effective scattering rate in a 
form similar to that of Eq (|A6|) : 

v n— p,7,... 

(5) 

where f n is the oscillator strength, hv n is the energy, and 
<t>n{x n ) is the Voigt functiorQ for the Lyman-n transition 
(i.e. Is — > np), where x n = (y — u n )/Avo. The first term of 
the series corresponds to Eq (|A6[> . The result of summing 
the first 10 Lyman series transitions is shown as the green 
dot-dashed line in Figure [2] Including these higher Lyman 
series photons only boosts the effective scattering rates by 
a factor of ~ 1.5 in the region that was optically thin to 
Ly/3 photons, which is due to the fact that the oscillator 
strength of np — > Is transitions -and hence the absorption 
cross section- decreases quite rapidly with n (e.g. for lOp — > 
Is, /io = 0.0016, which is ~ 0.004/ a ). Between 10~ 2 - 10 _1 
kpc from the edge of the HII region, the contribution from 
higher Lyman series photons is the largest and boosts the 
overall scattering rate by a factor ~ 15. This larger boost 
is due the fact that at these distances the gas is not yet 
optically thick in these higher Lyman series lines. Note that 
we have accounted for this boost in our calculation of tfs 
(which is shown Figure |3). 

(III) : We have so far assumed that the edge of the HII 
region is infinitely sharp. In reality, the thickness of the HII 
region is determined by the mean free path of the ionizing 
photons emitted by the quasar. Especially the X-Rays emit- 
ted by the quasar can penetrate deeply into the neutral IGM 
and produce a low level of ionization. However, the presence 
of a partially neutral region does not affect our results at 
all: although the distance from the edge of the HII region 
at which the IGM becomes optically thick is increased by a 
factor of 1/xm (where xhi is the neutral fraction), tfs is 
reduced by the same factor. 

(IV) : The approximation AT b (u) = -Tfs(v)T C mb/0- + 
z) is valid only when T ox <C Tcmb, or when n; > n u . The 2s 
and 2p populations are given by n 2s = ni = nu^ p / {A2 S is) 
and n 2p = n u = nH^ a /{A 2p i a ), respectively. That is, rt; > 
n u when F^ > (A 2s \ a / A 2p i s )Y a ~ 10~ 8 r a . Our assumption 
ni > n u therefore breaks down in the regions of interest only 
when r a J>10 -4 s~ 4 , which is orders of magnitude larger 
than the actual Lya scattering rate. 

(V) : Our calculation has ignored the hyperfine split- 
ting of the 2s and 2p levels that was illustrated in Figure [T] 



4 The oscillator strengths can be found Chapter 10.5 in Rybicki & 
Lightman (1979). The Einstein coefficients, and therefore the ab- 
sorption cross-sections, can be derived from the oscillator strength 
using Eq 10.79 of Rybicki & Lightman 1979. 
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Hyperfine splitting results in 3 allowed 2s- 2p transitions at 
v = 9.852, 9.876 and 10.030 Ghz, and the fine-structure ab- 
sorption cross-section has maxima at these frequencies. The 
most prominent of these 3 maxima reaches 75% of the value 
that o ne obtains when hyp erfine-structure splitting is ig- 
nored (|Dennison et al.ll2005h . Including hyperfine-structure 
therefore lowers the maximum observable brightness tem- 
perature by a factor of 0.75. 

5.2 General Detectability of the Fine-structure 
line during the EoR 

We have shown that the optical depth in the 2s-2p fine- 
structure line around quasar HII regions to be negligible be- 
cause of the low effective Ly/3, 7, 8, ... scattering rates. Early 
during the EoR, the IGM was exposed to a background Ly- 
man series photons generated by early gen erations of galax- 
ies b efore it was significantly ionized (e.g. iFurlanetto et al.l 
2006). In principle, one might imagine looking for fine- 
structure absorption in the IGM during this early phase of 
the EoR. Unfortunately, the effective Ly/3 scattering rate 
during this phase of the EoR will be low for the same rea- 
sons it is low around bright quasars, as we considered pre- 
viously: it is well established that Lya scattering rate plays 
an important role in determining the excitation tempera- 
ture of the 21-cm hyp erfine transition via the Wouthuysen- 
Field effect (see e.g. iMadau et all 1 19971 ; IFurlanetto et al.l 
2006). However, each soft-UV photon that redshifts toward 
(and through) the Lya resonance frequency scatters on av- 
erage iVscat ~ tgp ~ 10 5 — 10 6 times. In sharp contrast, 
each higher Lyman series photon only excites one atom into 
the 2s-state. The effective Ly/3 (and total) scattering rate 
is therefore lower than the Lya scattering rate by a fac- 
tor of ~ TGP- Since the Lya scattering rate in the IGM 
during the EoR is typically P a ~ 10" 12 - 10~ n s" 1 (e.g. 
IMadau et ai]|l997l ). the effective Ly/3 scattering rate is of 
the order Pg >e s ~ 10 -18 — 10 -16 s _1 , which barely competes 
with the recombination rate. From this we conclude that 
wherever the Wouthuysen-Field effect operates in the low 
density neutral IGM, the optical depth in the fine-structure 
line is expected to be tfs < 10~ 10 . 



the EoR. We have found that properly considering radiative 
transfer effects reduces the signal significantly compared to 
optically thin estimates, which radically complicates the de- 
tectability of the 3-cm fine structure transition. 

The main reason for this negative result is the large 
opacity of the (partially) neutral IGM to Lyman series 
photons. In such a medium, Lya scattering proceeds com- 
pletely different than Ly/3, 7, 5,... scattering: Lya photons 
scatter tgp ~ 10 s — 10 6 times. However, higher Lyman se- 
ries photons have a finite probability for being converted 
into Balmer, Paschen, .., etc photons at each scattering 
event, which subsequently propagate to the observer unob- 
structed. These higher Lyman series photons therefore scat- 
ter on average only 5 — 8 times f i]3.2[) before being destroyed 
(|Pritchard fc Furlanettd I2006T 1. Furthermore, only one of 
these 5 — 8 scattering events - defined as an effective scatter- 
ing event- indirectly excites atoms into their 2s-state (e.g. 
Ly/3 scattering corresponds to a sequence Is — > 3p — > Is, 
and does not affect the level population of the 2s-state). For 
these reasons, we found that the effective Ly/3 scattering 
rate was calculated reasonably well if one simply suppresses 
the incoming Ly/3 flux by a factor of e _T , in which r is the 
(frequency dependent) optical depth of the IGM to Ly/3 ra- 
diation (see Fig[5|. The same applies to Ly7, 5, ... radiation. 

Because of the reduced effective scattering rates, we 
found a substantial lower column of neutral atoms in the 
2s-state. This, in combination with the fact that the profile 
of the fine-structure absorption cross section is dominated 
by the damping wings (§(U), results in tfs ~ 10" 11 — 10 -10 . 
This is 5-6 order s of magnitude sm aller than what was orig- 
inally found by ISethi et all (|2007f ). and we conclude that 
the 3-cm fine-structure absorption line from neutral inter- 
galactic gas surrounding high-redshift quasars is presently 
undetectable. 
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6 CONCLUSIONS 

It has recently been proposed that neutral intergalactic 
atomic hydrogen gas may be detected in absorption in its 
3-cm fine-structure line (2si/2 — * 2p 3 / 2 ) against the Cosmic 
Microwave Background (CMB) out to very high redshifts. 
In particular, bright quasars that are luminous in the rest- 
frame soft-UV (10.2 eV < E 1 < 13.6 eV) may indirectly 
photoexcite the 2s state of neutral hydrogen gas in the in- 
tergalactic medium (IGM) during the Epoch of Reionization 
(EoR) to levels such that t fs ~ 10~ 5 for CM B photons pass- 
ing through these regions (|Sethi et alj|2007h . The resulting 
brightness temperature of ATj, ~ tens of mK could b e de- 
tected with existing radio telescopes (e.g. Carilli 2008). 

Motivated by this proposal, and by observational ef- 
forts to detect this transition, we have performed a detailed 
analysis of the transfer of Ly/3, 7, 5,... radiation, and have 
re-analyzed the detectability of the fine-structure line in 
neutral intergalactic gas surrounding bright quasars during 
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APPENDIX A: ANALYTIC CALCULATIONS 
OF LY/3 SCATTERING RATES 

Al The Optically Thin Limit 

In the regime where the gas is optically thin to Ly/3 radia- 
tion, the Ly/3 scattering rate is given by 



: 4-7T 



hv 



(Al) 



where Jiy) is the mean intensity in the UV radiation field 
near the Ly/3 resonance (in erg s _1 cm -2 Hz -1 sr _1 ), and 
opiy) is the Ly/3 absorption cross-section. The mean inten- 
sity J(v) relates to the luminosity of the central source (in 
erg s _1 Hz -1 ) as J{v) = 16 ^i r i , where r is the distance to 
the central source. Assuming that J(u) is constant across 
the Ly/3 resonance, which is extremely narrow in frequency, 
Eq [Al] simplifies to 



hv 



0, thin 



47rr 2 hvp 



J dv ap(v), (A2) 



where the photon energy hv, that was originally in the inte- 
grand, was taken outside of the integral since it barely varies 
over the limited range of frequencies where the integrand is 
not negligibly small. Using that J dv cp(v) = fp- 2 ^—, where 
fp = 0.079 is the Ly/3 oscillator strength, c is the speed 
of light, e and m e are the electron charge and mass, respec- 
tively (see e.g. lRvbicki fc Lightmanlll979l . their Eq 3.66), the 
total Ly/3 scattering rate simplifies further to 



P,thh 



h v fp ne 
inr 2 hvp m e 



(A3) 



The only quantity that remains to be determined is L„. Un- 
der the assumption that h v does not vary with frequency 
near the Ly/3 resonance frequency, h v = h/Av, where h is 



the total luminosity (in erg s x ) of the central source be- 
tween frequency i/ min and v max , and Av = (v max — Vm in ). 
S07 assume thalQ Np = 10 58 photon s _1 and Av ~ Q.lvp, 
and obtain 

'5 Mpc\ 2 



iWr) = 2.4 x 10"" s- 1 x(^)' 



(A4) 



where we used that h = Np x hvp = 1.9 x 10 47 erg s" 1 , and 
r = 5 Mpc (§ 13. ip . This scattering rate was plotted as the 
dotted line in Figure [2] 



A2 Including Radiative Transfer: Ignore 
Scattering 

When the assumption of working in the optically thin regime 
is dropped, the scattering rate (Eq IA1|) modifies to 



Pp (r) = 47T 



dv 3 -^Lop(y\\ - Hr/c}) e~ T ^ r) , (A5) 



where r(v,r) is the optical depth to photons emitted at fre- 
quency v to a, hydrogen atom at a distance r from the central 
source. Furthermore, the factor [1 — Hr/c] accounts for the 
fact that the Hubble expansion redshifts the photons by a 
factor of An = Hr over a distance r. Eq |A5l does not account 
for the fact that individual Ly/3 photons are scattered multi- 
ple times, which modifies the actual value of mean intensity 
(see §EHJ. 

For illustration purposes it is useful to rewrite Eq IA5I 
into a form similar to Eq IA1I 



Pg(r) 



L v fp ne 
Anr 2 hvp m e c 



x Sir), 



(A6) 



where S(r) (the 'suppression factor') is given by 



Sir) = fdx cj,(x)e- T(x - r) . (A7) 
V 71 " J 

Here, x denotes a dimensionless frequency variable of the 
form x = (y — vp)/AvD, where Avd = vthVp/c and 
Wth = *j2kT jm v . Furthermore, <j>{x) is the Voigt function 
(see Eq 10.77 of Rybicki & Lightman 1979. Note that this 
function obeys -\= J dx 4>{x) = 1. Hence, we recover Eq IA3I 
when we set r = 0.). The function 5(r) quantifies the inac- 
curacy in assuming that the absorbing gas is optically thin. 

We again assume that the edge of the HII region is 
infinitely sharp (as in § |4| and extends out to r = i?mi = 5 
Mpc (see § 15.11 roman numeral III). The optical depth r(r, v) 
is then 



r{v,r) 



ds nm(s)<rp(v[l — Hs/c]). 



(A8) 



Run 



We plot the suppression factor Sir) in Figure IA1I which 
shows that S(r) < 10" 2 for r > 10" 5 Mpc. The reason 
that Sir) drops so rapidly can be understood as follows: 

5 This assumption corresponds to a flux density that is 6 x 10 32 
which corresponds to an absolute AB magnitude 



erg s 



Hz 



of Mab = —30.45. For compari son, t he brightest z = 6 quasars 
have M1450 ~ -27.8 llFan et al.ll2006h . Given that quasars have 
a relatively fiat spectrum between A = Xp and A = 1450 A (rest- 
frame), we conclude that the quasar luminosity assumed by S07 
is likely too high by about a factor of ~ 10. 
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Distance from edge (Mpc) 

Figure Al. This figure shows suppression factor S(r), which 
quantifies the reduction of the Ly/3 scattering rate if one assumes 
an e _r suppression of the Ly/3 mean intensity. Because of the large 
Ly/3 absorption cross section the suppression is strong: S(r) < 
10 -2 for r > 10 -5 Mpc (see text). 



regardless of their frequency, photons traverse by definition 
an average optical depth r = 1 into the neutral intergalactic 
medium before being absorbed. For Ly/3 photons, this im- 
plies that the vast majority of photons are absorbed well be- 
fore they have redshifted into resonance. That is, t(x, r) ~ 1 
when x 3> 1 and therefore <fr(x) <fC 1, in the notation of 
Eq IA7I Only those photons that have redshifted into res- 
onance at the very edge can scatter near resonance (i.e. 
4>(x) = 1 while r(x,r) <C 1, also see £|3.2[) . This special 
region is very thin: the line-center optical-depth of neu- 
tral intergalactic gas in the Ly/3 line is to = (Tpfinmr = 
2.4 x 10 3 (T gas /10 4 )" 1/2 (r/kpc). That is, the line center op- 
tical depth approaches unity at r ~ 10 -3 kpc=10 -6 Mpc. 
The scattering rate that is obtained by combining Eq |A6llA"8l 
plotted as the dashed line in Figure [2] 



